Neel order in the two-dimensional S = ^-Heisenberg Model 

Utc Low 1 

1 Theoretische Physik, Universitat zu Koln, Ziilpicher Str.77, 50937 Koln, Germany 

(Dated: February 1, 2008) 

The existence of Neel order in the S = \ Heisenberg model on the square lattice at T = is 
shown using inequalities set up by Kennedy, Lieb and Shastry in combination with high precision 
Quantum Monte Carlo data. 



The ground state order of quantum spin systems, in 
particular the issue whether the ground state shows long 
range magnetic order, has attracted long and continuous 
interest. For the prototype of spin models, the antifcrro- 
magnetic Heisenberg model, the existence of Neel order 
at low temperatures was proved in the seminal paper of 
Dyson, Lieb and Simon [l| in 1978 for spin S > 1 and 
spatial dimension d > 3 and also for S = | and d > 3. 

Ten years later Kennedy, Lieb and Shastry Q showed 
that also for S = | and d = 3 Neel order in the ground 
state exists. 

The situation in two dimensions is different and more 
subtle, since the Mermin-Wagner-Hohenberg theorem 
forbids Neel order at finite T, leaving open however the 
possibility of Neel order in the ground state. The exis- 
tence of Neel order for the two-dimensional model and 
S > 1 was shown in [H, Q and later in Q by an indepen- 
dent derivation of the relevant inequality at T = 0. 

However the inequalities sufficient to show Neel order 
for S = 1 in the two-dimensional case are not sufficient 
to construct an analogous proof for S = g. Thus the 
case of S = ^ remains an open problem. Still it is pos- 
sible to derive inequalities concerning spin-spin correla- 
tions at short distances [H which are violated if Neel order 
is present. That is, with a minimum of numerical infor- 
mation, the question of Neel order in the ground state 
can be decided. 

The issue of this paper is to evaluate the spin-spin cor- 
relations of the two-dimensional S = | antiferromagnctic 
Heisenberg model at short distances and demonstrate 
that these results combined with the analytic expres- 
sions of Q show the existence of Neel order in the two- 
dimensional S = g antiferromagnetic Heisenberg model 
at T = 0. Such a study has become possible, due to the 
developement of high precision Monte-Carlo techniques 
over the last decade. 

In Ref.[|J Kennedy, Lieb and Shastry used data of 
Gross, Sanchez- Velasco and Siggia [|[ for a comparison, 
however these data clearly deviate from the results pre- 
sented here. The authors of [j| used a Quantum Monte 
Carlo method without loop updates and with discrete 
Trotter time (see below). Their data served only as a 
crude comparison to extrapolated Lanczos data and data 
produced by the Ncumann-Ulam method, which were 
the best algorithms to study the properties of the two- 
dimensional Heisenberg model in 1988. Today modern 
loop algortihms by far outreach both methods. 

As will be shown in the following an accurate evalua- 



tion of correlation functions at short distances is possible 
with modern Quantum Monte Carlo methods, which al- 
low us to compute expectation values at very low temper- 
atures and even though the short distance results have a 
certain finite size and finite temperature correction, these 
uncertainties are well controlled and allow to draw defi- 
nite conclusions. 

The approach and intention of this paper is diffrent 
from a completely numerical evaluation of e.g. the corre- 
lation length, which involves a calculation of correlations 
at long distances and an appropriate extrapolation to in- 
finite distances, which cannot be used as a proof of long 
range order in any rigorous sense. 

At first sight a "Quantum Monte Carlo algorithm" 
seems a puzzling concept, since an important step in 
any Monte-Carlo-method is the evaluation of Boltzmann 
weights for given energies of the system. For quantum 
models these energies are hard if not impossible to calcu- 
late. A key idea to make Monte Carlo methods applicable 
to quantum systems is to map the quantum model onto a 
classical model by introducing an extra dimension, usu- 
ally referred to as Trotter-time Q. 

In the first generation of algorithms this mapping 
was straightforwardly applied to the quantum Heisen- 
berg model. Though this allowed for a wealth of new 
studies of the finite temperature properties in one and in 
particular in two-dimensional systems, these algorithms 
had two major drawbacks, which became most evident 
at low temperatures. Firstly the extra Trotter dimension 
was discrctized, introducing the number of time slices as 
a parameter which had to be eliminated from the final 
results by an extrapolation. Secondly the update pro- 
cedure, i.e. the construction of new independent config- 
urations, was done locally. As a consequence one had 
to move through the lattice site by site several times to 
obtain a configuration independent of the starting config- 
uration and useful for a new evaluation of an observable. 

A first improvement was introduced by the so called 
loop-algorithms 0, which uses nonlocal updates similar 
to the Swcndsen-Wang algorithm for classical models. A 
second and important step towards high precision Quan- 
tum Monte Carlo techniques were algorithms which work 
directly in the Euclidian time continuum [8[ and require 
no extrapolation in Trotter time. For the algorithm Q 
used for the analysis presented here no approximations 
enter, and statistical errors are the only source of inac- 
curacy. 

Since this work intends to produce highly accurate 
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data it seems appropriate to assess the precision of the 
method by a comparison with exact results. The best 
candidate for such a comparison are the correlations of 
one-dimensional systems evaluated by the Bethe-ansatz 
with almost arbitrary precision up to distance seven [Io| 
and with results for finite chains from Ref. [HI]. This is 
done in the Appendix for chains of 400 sites at T=0.005. 

After these introductory remarks we now return to our 
actual goal, which is the two-dimensional system. Our 
starting point is a S = | Heisenberg model 



h = S X Sy 



(1) 
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with nearest neighbour interaction on a finite square lat- 
tice A with an even number of sites in every direction 
and periodic boundary conditions. 

The Fourier transform of the spin-spin correlation 
function at T = is given by 
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where 
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For the corresponding finite temperature expectation 
value of g q an upper bound f q was derived in [l L The 
T = limit of this bound was obtained in Ref. 4| and 
a direct proof of the bound at T = was given in Q. 
Following the notation and arguments of Kennedy, Licb 
and Shastry Q the inequality for d = 2 reads 



g q < f q tor q^Q 



(4) 



where /, 
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, E q = 2 - cos qi - cosg 2 , Q = 



(tt, 7r) and — eo is the ground state energy per site of the 
Heisenberg model Eq[T]on the lattice A. 

The fundamental idea is, that the existence of Neel 
order in the limit of infinite system size corresponds to a 
delta-function in the Fourier transform of the spin-spin 
correlation g q at Q. This means, if Eq. [¥] is integrated 
over the whole Brillouin zone one finds in the case of Neel 
order 



m 2 + 



J d 2 q f q > J d 2 q g q = S(S 



l)/3 
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where m 2 is the coefficient of the delta- function at Q. 

If there is no Neel order m 2 is zero. By numerically 
evaluating the integral over f q , and by using exact varia- 
tional upper and lower bounds on the ground state energy 
— eo one sees, that the above inequality and its analogon 
for d > 3 cannot be fulfilled with m 2 = and S > 1, 
which proves Neel order. 

Inequalities of type Eq. [5] are not sufficient to prove 
the existence of a nonzero m 2 for d = 2, 3 and S = 



but a new relation is obtained by multiplying g q by cos qi 
and again integrating over the Brillouin zone: 



d d qg q cosg 4 = (S* 2 5| i ) = -e0/3d 



(6) 



with i=l,2 for d=2 and i=l,2,3 for d=3 and Si the unit 
vector in i-direction and the value of the ground state 
energy form Ref. [Hi is e = 0.669437(5). 

Carrying out an analogous integral over f q and using 
again Eqf4]onc finds: 
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were the / + means the positive part of a function, which 
equals f, when f is positive and is zero otherwise. 

Again Eq[71 which is valid if no Neel order exists, was 
shown to be violated for d = 3 and S = h in Ref. [4] by 
using bounds on eo and thus the existence of Neel order 
was proved also for d = 3 and S = \ . 

For S = 7} and d = 2 one cannot construct a contra- 
diction by using only the ground state energy. Here more 
input from numerical data is needed. This can be incor- 
porated by multiplying g q by cos(mgi) with m = 2,3... 
and again integrating over the whole Brillouin zone: 



d qg q cos (mtft) = (S S mSz ) 



(8) 



with 1=1,2. 

Next, defining g(n) as 



1 - 
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and using again inequality 2] one constructs the follow- 
ing relations involving the correlation functions: 

f 1 - 

g(n) = / d2( l 2n + 2 2J ( -1 ) ro { cos ( m?1 ) + cos(mg 2 )} g< 
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Whenever the inequality Eq. [T0]is violated for a cer- 
tain n, a nonzero m 2 multiplying the delta-function at 
Q is needed and therefore the existence of Neel order is 
proved. 

The g(n) as defined in Eq[9] were calculated by the 
Quantum Monte Carlo method [3]. The results, displayed 
in table U show that the g(n) calculated by Quantum 
Monte Carlo cross the bound obtained by integrating 
over f q at n = 8. This is also depicted in Fig. [TJ Thus 
inequality Eq llOl is violated and Neel order must exists in 
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TABLE I: Bound obtained by integrating numerically over 
the right hand side of Eq 1 101 compared with g(n) for a 40 x 40 
lattice and different temperatures. 



TABLE II: Bound obtained by integrating numerically over 
the right hand side of Ea llOl compared with g(n) extrapolated 
for N=40,36,32,24 at T = 0.025. 



the two-dimensional antifcrromagnctic Hciscnbcrg model 
with S = \ at T = 0. 

There are three type of corrections to the data of table 
HI which need to be taken into account, but which, as 
we shall show in the following, do not change the above 
conclusion of a crossing of the curves at n = 8: 

(i) effects of finite temperature, 

(ii) effects of the finiteness of the system, 

(iii) statistical errors. 

In the following we comment on how these corrections 
modify the data. 

(i) The Quantum Monte Carlo data presented are at 
T > 0.005. The overall effect of finite temperature is to 
lower the absolute value of the correlations and therefore 
also the value of the g(n) . The effect of finite temperature 
is to shift the crossing of the bound and g(n) to larger n, 
or eventually to destroy a crossing completely. 

The functional dependence of the internal energy 
U(T), which up to an overall factor 3z (z = 2 is the 
coordination number of the two-dimensional square lat- 
tice) equals the correlation-function at distance o ne, has 
been determined for low T by spin wave theory [131 llij 
as 

U(T) = -e + bT 3 . (11) 

The coefficient is given in [H as b = ~ 0.2853626, 

so the correction for distance one is f» |l0~ 7 , which is 
two orders of magnitude smaller than the statistical error, 
(see point (iii)). 

For distances larger than one, we fitted the data as 
a function of temperature (taking the exponent of T as 



fit parameter) for T = 0.005, 0.025, 0.05, 0.075 and found 
the corrections due to finite temperature all of the order 
of 10~ 5 , which is the order of the statistical error. There- 
fore we do not give any finite temperature corrections. 

(ii) The absolute value of the correlations in the ther- 
modynamic limit are smaller than in systems of finite 
size. This means that the effect of finite system size is 
opposite to the effect of temperature. The finite size be- 
haviour of the ground state energy is well studied for 
the Hciscnbcrg model on the square lattice. Arguments 
originating from the quantum nonlinear sigma model de- 
scription [16| of the Heisenberg model to lowest order in 
system size give 

- eo = -eo(N) + , with c> (12) 

where —eo(N) is the ground state energy of a system of 
size N x N. Though the corrections are not substantial, 
they do effect the results, and taking into account, that 
the finite size errors in contrast to the finite temperature 
effects, might falsely lead to a crossing, we extrapolated 
the data for N = 24. ..40 using the functional dependence 
EqfT2l which we found well satisfied also for larger dis- 
tances. The results are shown in table HU One sees that 
the numeric values are changed but the crossing point is 
still at n = 8. 

(iii) We compute Ax = ^^ MC \/ i x2 ) — (x) 2 (where the 
observable x stands for the value of a correlation at a 
given distance, temperature and system size and Nmc is 
the number of Monte Carlo iterations), which is a reliable 
estimate for the statistical error of the mean value (x), 
since for the algorithm of Ref . [i[ the autocorrelation time 
is of order one and the Monte Carlo configurations are al- 
most independent. To assess the quality of our error anal- 
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FIG. 1: Bound on g(n) obtained from Eq|9]and g(n) for 24 x 
24, 36 x 36 and 40 x 40 at T = 0.025. 
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TABLE III: Correlations for a chain with N = 400 sites at 
T/J=0.005 compared with results from Ref. fTot | . 

ysis we also returned to the case of the one-dimensional 
antifcrromagnctic Hcisenbcrg model ( see Appendix ) and 



compared results with independent streams of random 
numbers. 

To calculate an upper limit to the errors of g{n), the er- 
rors of the correlations where added up ( being evaluated 
with the same configurations, they are not independent). 

To conclude, the error analysis shows that the short 
range correlations entering Eq[9] were determined with 
sufficiently high accuracy to prove the existence of a 
crossing of the bound and the Quantum Monte Carlo 
data for g(n) at n = 8 and therefore to show the 
existence of long range order. 

Appendix 

(1) In this Appendix we list the correlations of a one- 
dimensional Heisenberg model with periodic boundary 
conditions and chain length N = 400 at T=0.005 com- 
pared with results of Ref. [H| for infinite chain length and 
T = 0. 

For the internal energy U (T) of the Heisenberg chain 
the temperature dependence for low T is U(T) = —e\ + 
aT 2 with the ground-state energy = 0.4431471804 
for 400 sites and ej = — \ + In 2 for the infinite size 
system fl7j . and the coefficient a = -| given in Ref. 
0, Ell- This means that the correction for the corre- 
lations in tabl dllll due to finite temperatures are of the 
order of 10~ 5 . 

(2) The exact values of the correlation functions 

[111 . l20l ] for distance one and two at T = for a chain 
with 400 sites are {S$Sf) 40 o = -0.147717441765735 and 
(S$Sl) 4 oo = 0.0606813790491800. The above data show 
that the error analysis concerning statistical errors and 
finite temperature effects is consistent. 
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